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PREDICTION OF FLAME VELOCITIES OF HYDROCARBON FLAMES 
By Gordon L. Dugger and Dorothy M. Simon 


SUMMARY 

The laminar -flame -velocity data previously reported hy the Lewis 
laboratory are surveyed with respect to the correspondence between 
experimental flame velocities and values predicted by semitheoretical 
and empirical methods. The combustible mixture ’ variables covered are 
hydrocarbon structure (56 hydrocarbons), equivalence ratio of fuel-air 
mixture, mole fraction of oxygen in the primary oxygen-nitrogen mixture 
(0.17 to 0.50), and initial mixture temperature (200° to 615° K) . The 
semitheoretical methods of prediction considered are based on three 
approximate theoretical equations for flame velocity: the Semenov 

equation (bimolecular) , which is primarily based on conductive heat 
transfer between the flame and the reactants; the Tanford-Pease equa- 
tion, which is based on the diffusion of chain carriers of the oxidation 
reaction into the reactants; and the Manson equation, which is a modi- 
fication of the momentum-pressure -drop equation that does not include 
chemical kinetics. In each equation a semiempirical factor is used to 
bring the predicted values for a given variable and fuel into the best 
average agreement with the data, so that the variation in the relative 
prediction of the individual datum points may be considered. For the 
resulting semitheoretical equations, it is assumed that thermal equilib- 
rium is attained at the end of the flame zone, and values for the trans- 
port properties are estimated by extrapolation and simple additive rela- 
tions. The empirical relations between combustible mixture variables 
and flame velocity are based on the usual methods of correlation. 

Within these limitations, the results may be generalized as follows 

1. The three semitheoretical equations predict relative flame 
velocities reasonably well, generally with mean deviations of 2 to 
15 percent. 

2. Considering only the bimolecular fuel-oxygen reaction and with 
low-temperature activation energies, the Semenov equation may be used to 
give good relative predictions when a "steric factor" is determined 
semiempirically . 

3. The average- " specif ic rate constants" (k a ) obtained with the 
Tanford-Pease equation for molecular fuel— active particle reactions 
behave as Arrhenius rate constants only with the data on the effect of 
initial temperature, not with the equivalence-ratio or. oxygen- 
concentration data. The relative predictions obtained by this equation 
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are not very sensitive either to the temperature dependence assigned to 
'the diffusion coefficients or to the recombination factor computed for 
hydrogen atoms. Regarding active particles considered, better results 
are generally obtained for hydrocarbon flames when the hydroxyl radical 
and the hydrogen and oxygen atoms (OH, H, and 0) are considered rather 
than H alone; this is particularly true for studies of the effect of 
equivalence ratio. 

4. For the equivalence-ratio and oxygen-concentration data, con- 
siderably better relative predictions are obtained from the Manson 
equation when the pressure drop across the flame front is considered to 
be due to H, OH, and 0 than when only H atoms are taken into account. 

5. The example presented by the ethylene data shows that an empiri- 
cally determined rate constant or proportionality factor from equivalence 
ratio data at atmospheric pressure and room temperature may predict the 
effect of initial temperature or oxygen concentration within approxi- 
mately 20 percent by the Semenov or Tanford-Pease equation and within 
approximately 30 percent by the Manson equation. 

6. For engineering applications, the effects of the parameters 
studied could be estimated just as satisfactorily, and more easily, by 
one or another of the empirical correlations indicated, as compared with 
the three semi theoretical equations considered. However, the use of the 
semi theoretical equations in some cases reduces the number of constants 
required. 


INTRODUCTION 

The ability to predict flame velocities of fuels is of growing 
importance in the field of aircraft propulsion, since a correlation has 
been found between combustion efficiency of a ram-jet burner and the 
laminar flame velocity of the fuel (reference 1). The prediction of 
flame velocities is difficult for three reasons, (l) There is no com- 
plete, rigorous theory which can be readily applied. There are, however, 
a number of approximate equations in the literature which approach the 
problem of flame .propagation from various viewpoints. (2) There are no 
data on the kinetics of the oxidation process under flame conditions - 
and very few data on transport properties at high temperatures. (3) 
Different methods of flame -velocity measurement give different values, 
so that it is difficult to compare data from different sources. The 
uncertainty in measurements made by a given method is of the order of 
5 percent. 

In this report flame -velocity measurements made at the NACA (ref- 
erences 2 to 7) for different hydrocarbons, initial temperatures, and 
compositions are used with semi theoretical and empirical methods of- 
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flame-velocity prediction to show the correspondence between the meas- 
ured velocities and the predicted velocities. The semi theoretical 
methods are based on the Semenov equation (reference 8), the Tanford- 
Pease square- root law (reference 9), and the Manson equation (refer- 
ence ll) . These three equations were derived from different models of 
the flame -propagation process. The Semenov model is essentially a 
thermal model which includes chemical reaction kinetics; the Tanford- 
Pease model is based on the diffusion of chain carriers of the oxidation 
reaction; and Manson used a modification of the momentum-pressure-drop 
equation which does not include chemical kinetics. The empirical rela- 
tions for the effect on maximum. flame speed of hydrocarbon structure, of 
initial mixture temperature, and of oxygen concentration were based on 
the usual methods of correlation. 

It is recognized that there are other approximate theoretical equa- 
tions which may give as good or better predictions of the NACA experi- 
mental results. These particular equations were chosen because they 
exemplify three different approaches to the problem and because the 
Semenov and Tanford-Pease equations were used in previous NACA papers 
and many of the calculations had already been made individually. It is 
also recognized that there are many flame -speed data in the literature 
which could be used in such studies; the present paper is confined to 
NACA data because many of the calculations had been made, and because it 
was desirable to avoid the complications arising from comparing data 
obtained by many different methods and techniques. 


NOMENCLATURE 

The following nomenclature is used in this paper: 
a fuel concentration, molecules of fuel per cm 3 of mixture 
b oxygen concentration, molecules of oxygen per cm 3 of mixture 

Bjl term near unity arising from radical recombination 

Cp molar heat capacity at constant pressure, cal/mole-°K 
Cp specific heat, cal/g-°K 

Cp mean specific heat, Tq to Tp, cal/g-°K 

D diffusion coefficient, cm^/sec 

diffusion coefficient of i x active species into unburned gas at 
initial mixture temperature, cm^/sec 
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activation energy, kcal/g-mole 

fitted activation energy for a group of data points, kcal/g-mole 

empirical proportionality constant between experimental flame 
velocity and value predicted by the Manson equation 


empirical K from Manson equation when Ap « l/2 pg 



empirical K from Manson equation when 

a ,lJ -l. d oh ^ D o\ T o 
Ap - l/2 ^p H + p QH + P 0 d h J T f 


average Kg or for a group of data points 

specific rate constant, cm^/molecule-sec 

/ 

weighted mean kj_ for the three active particles H, OH, and 
0, each reacting with fuel molecules, cm^/molecule-sec 

average k a for a group of data points, em^/molecule-sec 

specific rate’ constant for reaction between fuel molecule and 
ith active particle, cm^/molecule-sec 

total concentration of gas at mean combustion- zone temperature, 
molecule/cm^ 

molecular weight 

molecularity of flame reaction 

J-l* 

mole fraction of i • component 

total number of molecules of HgO and COg in products per molecule 
of fuel by stoichiometric relation 

moles of reactants per mole of products from stoichiometric 
equation 

steric factor, or probability factor, from the expression 
k = PZe-^ 1 

average P for a group of data points 
total pressure of mixture 
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p i mole fraction of the active particle in burned gas 
Q mole fraction of potential combustion product in unburned gas 

Q* mole fraction of fuel in unburned gas 

R universal gas constant, kcal/(g-mole) (°K) or ergs/(g-mole) (°K) 

T absolute temperature, °K 

computed equilibrium flange temperature, °K 

U flame velocity, cm/sec; implies maximum flame velocity with 
respect to equivalence ratio, except in discussion of the 
variation of equivalence ratio 

w reaction rate, molecules reacting per cm^/sec 

Z collision number; the number of molecular collisions per second 
when the concentration is one molecule of each type per cm^ 

a mole fraction of oxygen in oxygen-nitrogen portion of mixture, 

o 2 /(o 2 +n 2 ) 

& percentage mean deviation in the ratio of predicted flame velocity 
to the experimental value for a given group of data 

i] viscosity of mixture, poise 

~ viscosity of i ttl component, poise 

6 m ratio of mean reaction- zone temperature to initial temperature 

X thermal conductivity, cal/ ( cm 2 ) ( sec ) ( °K/ cm) 

p density of mixture, g/cm 3 

a collision diameter, cm 

cp equivalence ratio, fraction of stoichiometric fuel-oxygen ratio 

Subscripts : 

eff effective value at mean reaction temperature 
f at flame temperature 
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m at mean reaction-zone temperature 

0 initial conditions 


APPROXIMATE THEORETICAL EQUATIONS USED 


The Semenov Bimolecular Equation 


Zeldovich and Frank -Kamenetsky obtained an approximate solution 
for the rate of flame propagation from the differential equations for 
heat conduction and fuel concentration change across the flame front. 
This approximate solution was presented in detail by Semenov (refer- 
ence 8). The general approximate solution maybe written (see refer- 
ence 8, pp. 31, 48, and 49) : 


U = 


1 



a o p o c P 


w dT 



a) 


where, for a bimolecular reaction between unlike molecules 


T-f 


w dT - 


k(a) (b) dT 


. . P ~E/RT f 


- a e ff b e ff PZ 


E 


( 2 ) 


For lean mixtures, cp<l, a.Qff and b e jj are computed by (fol- 
lowing the assumptions of reference 8, p. 45) : 


RT* 


-Lq -g 

a eff =a O^ T f -T 0 


b eff sb 0f 


f L 


RT* 


1 - cp\l - 


E 


T f T 0, 


(3) 


For rich mixtures, cp > 1 : 


x eff =a 0 t 


£ 

'f L 


RT, 




E 


W j 


NACA RM E52J13 


7 


RT f 2 

b eff = b 0 T^T Tf -T 0 ^ 

For stoichiometric mixtures, equations (3) and (4) are identical. 
All of the factors in equations (l) through (4) can, with the exceptions 
of P and E, be estimated by the extrapolation of thermodynamic tables 
(references 12 and 13) or calculated by the relations 

T = I!) \ N i 

i 

(Values calculated for combustion-product mixtures by this simple, addi- 
tive method were within 1 percent of values calculated by the method of 
reference 14. ) 


C P - sp Cp,i N ± 




A = (Cp + j R) tj/M (reference 13) 

D s 1.336 T|/ p (reference 15) 
p = p/MRT 

\ 2 / L ^ 


> (5) 


J 


Flame -velocity predictions by the Semenov equation are evaluated 
herein as follows. Activation energies from low-temperature reactions 
and measured flame velocities are substituted in equation (l) and a 
steric factor P is calculated for each experimental flame velocity. 
These P's are then averaged to give P for the group of data points. 
The ratio of the predicted flame velocity to the experimental velocity 
is calculated as (P/P)l/2. The average deviation of the ratios from 1 
is considered to be a measure of the accuracy of the predicted flame 
velocities. (The word accuracy is used to denote the degree of corre- 
spondence between flame velocities calculated by the methods described 
and measured flame velocities.) 


The data for various hydrocarbons consisted of maximum flame 
velocities (maximum with respect to equivalence ratio) of 56 pure hydro- 
carbons at 298° K and atmospheric pressure (references 2 and 3) • These 
data were obtained by the NACA tube method. The hydrocarbons included 
straight and branched chain members of the alkane, alkene, alkadiene, 
alkyne, and cycloalkane series and benzene. No calculations are made by 
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the Semenov equation for these data, because the work required by this 
treatment would not be warranted in view of the slight differences in 
flame velocities observed for most of the hydrocarbons studied. However, 
it has been shown by more approximate calculations (reference 16) that 
relative predictions would be satisfactory. 

The initial-mixture -temperature data included maximum flame veloc- 
ities of methane-air, propane-air, and ethylene-air mixtures obtained by 
a Bunsen burner method (total area method, outside edge of shadow cast 
by f lam e cone, reference 5) at seven initial mixture temperatures -ranging 
from 200° to 617° K- The following activation energies, in kilocalories 
per gram-mole, are reported in references 17, 18, and 19, respectively: 
methane, 51; propane, 38; ethylene, 40. These values are used in the 
Semenov equation to evaluate the flame -velocity predictions. The accu- 
racy of the predicted flame’ velocities is 2 to 3 percent. Specific 
values axe tabulated in table I under "Semenov Equation" . 

Two kinds of mixture -composition variables are studied. The first 
is the change of hydrocarbon concentration in air over an equivalence - 
ratio r ang e of 0.7 to 1.3. ‘Flame-velocity values by the tube method for 
ethylene-air and pentane-air at 298° K (reference 4) , and Bunsen burner 
values (total area method, outer edge of cone shadow) for methane-air 
at 307° K and propane -air at 302° K (reference 5) are used. When the 
effect of changes in equivalence ratio on the flame velocity is computed 
by the Semenov equation, it is found that the approximations used in 
equations (3) and (4) are not very consistent for the region near 
stoichiometric, particularly for equivalence ratios between 0.95 and 
1.05. This fact is shown by the difference between the Semenov curves 
in figure l(a), where points in the range 0.95 to 1.05 ape included in 
plotting the curve, and figure l(b) , where these points are .omitted. 

If this critical region is omitted, the accuracy of the predicted flame 
velocities over the equivalence range is 3 to 7 percent. Specific 
deviations are listed in table I. > 

The second composition variable studied is the ratio of oxygen to 
oxygen plus nitrogen a = 0 2 /(0 2 +N 2 ). Maximum flame velocities for 

several a ratios were available for the following fuels and ranges of 
a: propane, 0.17 to 0.50; ethylene, 0-17 to 0-35; and 

2,2,4-trimethylpentane, 0.21 to 0.50. The accuracy of prediction of 
these data is 3 to 17 percent (table I) . 

It may be noted in table I that flame-velocity predictions by the 
Semenov equation for ethylene flames are accurate to 4 percent for 
variations in both composition and initial temperature. For these 
predictions (as discussed previously) , each type of data was considered 
separately, and steric factors were calculated from the data. Now 
usually, when it is desired to predict flame velocities, only one type 
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of data is available - for example, fuel-air ratio data at atmospheric 
pressure and room temperature. It is therefore important to know how 
well the effects of initial temperature and oxygen concentration can be 
predicted from such data. The ethylene data are used as an example. 

The steric factor P calculated from the equivalence -ratio data 
and the low -temperature activation energy of 40 kilocalories per gram- 
mole are used in the Semenov equation to predict flame velocities over 
the ranges of initial temperature and oxygen concentration covered by 
the experimental data. In both cases the predicted flame velocities 
deviated from the normalized measured flame velocities by an average of 
+14 percent, the maximum deviation being +24 percent for the oxygen- ' 
concentration data. (Because different methods of measuring flame 
velocity were used, differences in values for mixtures of the same com- 
position and temperature were found. For these calculations all flame 
velocities were normalized by a simple ratio factor to the values for 
the initial-temperature data.) 


Tanford-Pease Equation 

Tanford and Pease (references 9 and 20) equated the amount of 
product formed in the combustion zone by a second-order reaction between • 
fuel molecules and hydrogen atoms (or other active particles such as 
hydroxyl radicals or oxygen atoms) to the amount of product formed at 
the flame front by conversion of the fresh gas expressed in terms of 
initial conditions and flame velocity. An approximate solution for flame 
velocity from this equation was given: 

Three methods of evaluating the predictions of the Tanford and 
Pease equation are used herein. For all three, the following calcula- 
tions are the same: (a) flame temperature T f and active particle 

concentrations Pj_ are calculated assuming adiabatic thermal equilibrium 
by a matrix method (reference 21) ; (b) the mean combustion- zone tempera- 
ture is assumed to be 0.7 Tf (reference 9); -(c) diffusion coefficients 
D m are calculated from Dm “ ^i (0.7 T f/ T 0r* j where J)± .is the 
diffusion coefficient at intial temperature calculated by the Stefan- 
Maxwell equation (reference 22) ; (d) the recombination factor .is 

calculated by the method of Tanford (reference 10) for the hydrogen 
atoms and is assumed to be one for OH and 0; (e) the ratio l^Q'n/Q is 

calculated from a knowledge of the over-all oxidation process and the 
initial concentrations of reactants. 
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The three methods of evaluation differ in the calculation of kj_. 
For the first method only one chain carrier, the hydrogen atom, is con- 
sidered. For this case kg values are calculated from single-point 
flame-velocity determinations hy equation (6). For the second- method, 
H, 0, and OH are considered to he the chain carriers and 

c 


^i^i^m ~ , 

- k a 


Bj 


r PH D H 

V % 


B0H D 0H + Bo 11 © 


( 0.1 T4 


,1.67 


(7) 


Single-point k a values are calculated. For these two methods, the 
ratio of the predicted flame velocity to the experimental flame velocity 
is calculated as (TL^/k^) and the average deviation of the ratio 

from 1 is used as a measure of the accuracy of the predictions. 


For the third method, an Arrhenius type temperature dependence of 
the rate constant is used: 


ka - R a Z a ex P (- E a/ R (°‘ 7 T f ) ) 


( 8 ) 


The following method is used to calculate an activation energy E a for 
best fit from the initial-temperature data. Several values of E a are 
substituted in the equation and for each a value of (P/P) is cal- 

culated. Then 6, the percentage mean deviation in (p/p)V R fob" the 
group, is plotted against the E a used. The "fitted" value of E a is 
the one which gives the minimum &. All other calculations by the third 
method are made by equations (7) and (8), measured flame velocities, and 
E a from initial -temperature data. The accuracy is calculated as the 
mean deviation of (P/P)V R from 1. 


The predictions of the Tanford-Pease equation are evaluated by 
methods 1 and 2 for the maximum flame velocities of a series of hydro- 
carbons (reference 23) . Table I shows that the accuracy of prediction 
of flame velocity is 3.3 percent considering hydrogen atoms (kg) only 
and 3.5 percent considering all three active particles (ka). These 
mean deviations are small enough to indicate either that the "rate con- 
stants" are not temperature -dependent, in which case the steric factors 


are of the order of 10 , or that 



is nearly con- 


stant for the hydrocarbons studied, which is probable. If the activation 
energy is- of the order of 10 kilocalories per gram -mole, the steric fac- 
tor is of the order of 10 _E . 


All three methods of evaluation are used with the initial- 
temperature data. The accuracy of the flame velocity predictions is 
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9 to 12 percent (table I) for the first method (kg) and 6 to 8 percent 
for the second method. The use of method 3, which includes the tempera- 
ture dependence of the rate constant, improves the prediction so that 
the accuracy is about 2 percent. The calculated activation energies and 
steric factors are: for methane, 17 and 0.19; and for ethylene, 18 and 

0.21. The calculated activation energy for ethylene is high when com- 
, pared with the experimental activation energy for the reaction 
H + C2H4— >C2H5, which is 2.6 kilocalories (calculated from data in 
reference 24) . The relative predictions are not appreciably changed by 
assuming Bg » 1 or by assuming D m cc T m which was used in reference 7 
instead of DjjjCOq/p. 

If the entire equivalence -ratio ranges which were covered experi- 
mentally by the flame -velocity data are used for the prediction of 
velocities by the Tanford-Pease method, the accuracy of prediction is 
very poor as illustrated in figure 1(a) for ethylerie. If the equivalence 
range is limited to one near the maximum flame velocity (for example, 
cp a c p mnv ^Q. 15) , the accuracy of the predictions is 1 to 6 percent 

(table I) when all active particles are considered. However, if only 
hydrogen atoms are used the accuracy is very poor - 17 to 65 percent. 

The variations of flame velocity with oxygen concentration were 
calculated for all active particles. The accuracy of prediction is 4 to 
6 percent. 

The ethylene data are_again used as an example to show how well an 
average "rate constant" (1% from method 2) from the fuel-air ratio 
data predicts the effects of temperature and oxygen concentration on the 
maximum flame velocity. Predicted flame velocities for various tempera- 
tures differ from the normalized experimental values by -9 to +13 percent 
with an average deviation of 6 percent. Predicted flame velocities for 
various oxygen concentrations differ from experimental values by +8 to 
+22 percent with an average deviation of 12 percent. 


The Mans on Equation 

The momentum relationship between. the flame (or combustion wave) 
velocity and the pressure drop across a plane, steady-state flame front 
may be expressed in the form (reference 25, pp. 24lX: 



a p y / 2 

p o- p f/ 


O) 
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Manson (reference ll) suggested that the small pressure drop could he 
caused by the projection of hydrogen atoms into the uriburned gas. Because 
the H atoms would recombine to H2 at the uriburned gas temperature, 

Ap was assumed to be one-half the equilibrium H atom pressure reduced 
to uriburned gas temperature. In the present paper, Ap is calculated as 
1 T 0 1 ( Doh Do \ To 

both g p H and 2 yPH + P 0H pj" + p 0 % J T^" The second method of 

calculating Ap takes into account three of the lightest and most 
abundant species which might be considered active particles and assumes 
that all of these would recombine in pairs. The average proportionality 
factor for a group of data is calculated by dividing the experimental 
flame velocity by the right-hand side of equation (9) for each point and 
averaging the quotients. These average empirical proportionality factors 
are designated Kg for the first method and Kg for the second. The 
accuracy of the predicted flame velocities is considered to be the 
average deviation of %/% or Kg/Kg from 1. For the maximum flame 

velocities of the different hydrocarbons, the accuracy of prediction by 
the Manson equation for either method of evaluation is about 4 percent 
(table I) . For the first method (Kg) , the proportionality factor is 
0.5; that is, the predicted value is twice the experimental value. 

For the initial-temperature data the predicted flame velocities 
differ from the experimental velocities by 10 to 14 percent. When the 
total range of the equivalence -ratio data is 'considered, the accuracy of 
prediction of flame velocities by the Manson equation and either of the 
methods of evaluating Ap is poor . If the equivalence range is limited 
as for the Tanford -Pease method and with all active particles, predic- 
tions accurate to 2 to 4 percent are obtained (table l).> The accuracy 
of prediction of flame velocities for mixtures containing different 
concentrations of oxygen varies from hydrocarbon to hydrocarbon. For'’ 
2,2,4-trimethylpentane, the accuracy is about 1 percent; for ethylene, 

5 percent; and for propane, 15 percent. 

The ethylene data are again considered in order to determine the 
accuracy with which the Kg from the equivalence -ratio data predicts 
the effects of temperature and oxygen concentration on flame velocity. 

The deviations for the initial-temperature data varied from -9 to +45 per- 
cent with a mean deviationof -27 percent. For the oxygen data the 
deviations varied from +1 to +19 percent with a mean value of +10 percent. 


Empirical Correlations 

A correlation was obtained (reference 26) by which maximum flame 
velocities of various hydrocarbons in air could be predicted with an 
average deviation of 2 percent. The predicted flame velocity is calcu- 
lated from the sum of the contributions of various H-C bonds in the fuel 
molecule according to the following relationship: 
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U = n a k a + n b k b + + . . . 

where N A , N B , Nq, N b , N b , N b , Nq, and Ng- are the numbers of methane, 

primary, secondary, tertiary, alkene, alkyne, cyclohexyl and aromatic 
C-H bonds, respectively, per unit volume of hydrocarbon-air mixture, 
and K a , K b , K^, • . • are the flame-speed coefficients of these bonds. 

For the special cases of C-H bonds on carbon atoms placed alpha to the 
alkyne C = C bond, a factor M = 0.96 was introduced into terms repre- 
senting these alpha bonds. Correlation coefficients (table II) estab- 
lished from 34 hydrocarbons excluding ethylene then gave an average 
deviation of 1.9 percent in the ratio of predicted to measured flame 
velocity. 

A correlation was also' observed between the equivalence ratio for 
maximum flame velocity and the total bond dissociation energy of the 
fuel. As shown, in table III, the total bond dissociation energy per 
unit volume of the hydrocarbon-air mixture corresponding to the maximum 
flame velocity is nearly constant, with an average deviation of only 
0.9 percent from the average value. This deviation is less than one- 
third of that obtained by simply assuming that the maximum will occur 
at an average equivalence ratio of 1.15. Bond energy calculations for 
37 hydrocarbons similar to those presented in table III had an average 
deviation of 1.6 percent. 

For the equivalence -ratio data, linear relationships were found 
between flame velocity and the logarithm of the equivalence ratio, cp, 
where cp < c p mRy . An interesting variation is the plot of flame velocity 

against the total bond dissociation energy of the fuel per unit volume 
of mixture in figure 2 for pentane, ethylene, and propyne mixtures with 
air. It may be seen that the correlation is linear for the lean mix- 
tures and that the line extrapolates to the lean limit for flame propa- 
gation (where U = 0) as determined in the 1-inch flame tube. For these 
data, determined by the tube method for three hydrocarbons, good. esti- 
mates of the maximum flame velocities are obtained by reading from the 
curves at a total bond dissociation energy of 1.75 kilocalories per 
liter. This constant bond-energy value of 1.75 kilocalories per liter for 
this correlation does not equal the bond energy of the hydrocarbon con- 
centration corresponding to the maximum flame velocity which was pre- 
viously calculated to be 1.96 kilocalories per liter. Work with other 
flame -velocity data shows that both the position of the horizontal line 
and the slopes of the lines for specific hydrocarbons depend somewhat 
on the method of measurement of flame velocity. 

Empirical equations predicting the effect of initial mixture tem- 
perature on flame velocity with mean deviations of approximately 2 per- 
cent for the range from 200° to 615° K are (reference 5) : 
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for methane 


U « 10 + 7.40X10 T, 


-5 m 2.23 


for propane 


U - 10 + 3.42 X10 T, 


-4 m 2.00 


for ethylene 

U = 10 + 2.59 X10" 3 T 0 1,74 

Empirical relationships were found which predict the effect of 
oxygen concentration (for the experimental range covered) and limited 
temperature (311°-422° K) on the flame velocity for the three hydro- 
carbons studied (reference 6) : 

for 2,2,4-trimethylpentane 

U = 0.133 T 0 1 ' 40 (a-0.120) 


for propane 

U = 0.766 Tq 1 * 16 (a-0. 133) 

for ethylene 

U = 0.998 Tq 1 ’ 18 (a-0. 133) 

The accuracy of prediction was about 3 percent for ethylene and 
2,2,4-trimethylpentane and 6 percent for propane. 


CONCLUDING REMARKS 

The semi theoretical calculations of this paper are dependent on the 
assumption that thermal equilibrium is attained at the end of the com- 
bustion zone. Actual flame temperatures and product concentrations 
would differ from equilibrium values if (a) equilibrium is not attained, 

(b) any appreciable reaction takes place at the initial temperature, or 

(c) there is a large chain branching term affecting the radical distri- 
bution in and ahead of the flame zone. These calculations are also 
dependent on the inherent assumptions of the approximate theoretical 
equations and the methods used to calculate transport properties. Within 
these limitations, the results of the calculations may be generalized 

as follows : • 
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1. The Semenov bimolecular equation (thermal mechanism), the 
Tanford-Pease square-root law (active -particle diffusion mechanism), and 
the Manson modification of the momentum-pressure -drop equation (including 
diffusion of active particles, hut not chemical kinetics) will all predict 
relative changes in flame velocity caused by changes in hydrocarbon 
structure, initial temperature, equivalence ratio (limited range) , or 
oxygen concentration- The accuracies of the predictions are summarized 
in table I- 


2. Low-temperature activation energies may be used in the Semenov 
equation to give good relative predictions of flame velocities for var- 
iations in initial temperature, equivalence ratio, and oxygen concentra- 
tion- Only the bimolecular fuel-oxygen reaction is considered in this, 
paper; other combinations such as fuel and hydrogen atoms should be 
investigated to determine whether the kinetics indicated for such other 
reactions are not more probable than for the fuel-oxygen case- 

3- The average "specific rate constants" (kgX obtained with the 
Tanford-Pease equation for molecular fuel .- active particle reactions 
behave as Arrhenius rate constants only with the data on the effect of 
initial temperature, not with the equivalence -ratio or oxygen-concentra- 
tion data. The relative predictions obtained by this equation are not 
very sensitive either to the temperature dependence assigned to the 
diffusion coefficients or to the recombination factor computed for hydro- 
gen atoms; these factors do appreciably affect absolute predictions of 
flame velocity or, conversely, determinations of steric factors from 
experimental flame velocities- Regarding active particles considered, 
better results are generally obtained for hydrocarbon flames when H, 

OH, and 0 are considered rather than H alone; this is particularly true 
for studies of the effect of equivalence ratio. 

* 

4. For the equivalence -ratio and oxygen-concentration data, con- 
siderably better relative predictions are obtained from the Manson equa- 
tion when the pressure drop across the flame front is considered to be 
due to H, OH, and 0 than when only H atoms are taken into account. 

Further thought should be given as to how this pressure drop should 

be computed, taking into account any net diffusional flow of any 
component between the flame front and the unburned gas. 

5 . The example presented by the ethylene data shows that an empir- 
ically determined rate constant or proportionality factor from 
equivalence -ratio data at atmospheric pressure and room temperature may 
predict the effect of initial temperature or oxygen concentration within 
approximately 20 percent by the Semenov or Tanford-Pease equation and 
within approximately 30 percent by the Manson equation. 

6. For engineering applications, the effects of the parameters 
studied could be estimated just as satisfactorily, and more easily, by 
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one or another of the empirical correlations indicated, as compared 
with the three semitheoretical equations considered. However, the use 
of the semitheoretical equations in some cases reduces the number of 
constants required. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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TABLE I - PERCENTAGE MEAN DEVIATIONS IN RATIO OF PREDICTED TO MEASURED FLAME VELOCITY 


Type of data 

Empirical correlation 

Semenov 

bimolecular 

equation 

Tanford-Pease 

equation 

Mans on 
equation 

Form of correlation 

Deviation 

(percent) 

Using 

k H 

Using 

k a 

Using 
E from . 
Tq data 

Using 

% 

Using 

K s 

Hydrocarbon structure 8. 

u = n a k a + n b k b + . . . 

1.9 


3.3 

3.5 


3-6 

3.8 

Initial temperature^ 

(200° to 615° K) 

Methane 

Propane 

Ethylene 

U = a+bT 0 n 

2.1 

1.6 

1.4 

2.7 
2.0 

1.8 

11.8 

12.5 

9.4 

7.8 
6.1 

5.8 

0.8 

1.8 

10.3 

10.0 

12.0 

10.1 

13.0 

13.8 

Equivalence ratio 
(range limited, see text) 
Methane, 307° Kr 
Pentane 8, 

Propane, 302° K b 
Ethylene 8, 

U = a ( log B) + b 

4.9 

2.2 

7.4 

3.0 

4.0 

5.0 

25.8 

64.8 

16.6 

5.7 

3.8 
0.6 
3.7 

7.0 

3.3 

23.6 

39.9 

15.1 

10.0 

1.8 

3.6 

1.2 

1.8 

Oxygen concentration 0 
Propane, 311° K 
422° K 

Ethylene, 311° K 
422° K 

2,2,4-trimethyl- 311° K 
pentane 422° K 

|u a aTg 11 (a-b) 

Ju = aT 0 n (a-b) 

| U = a (a-b) 

U = aT 0 n (a-b) 

J 5.5 

J 3.2 

2.0 

2.7 

2.6 

17.3 

12.4 
4.0 
3-0 
6.4 
7.3 


6.4 
5.2 
4.9 
4.6 

4.5 
4.0 


16.5 

12.7 

7.5 

6.8 

7.3 

8.2 



Data determined by tube method at room temperature (reference 2) . 

Data determined by Bunsen burner method, shadow cone (reference 5) . 
c Data determined by Bunsen burner method, schlieren cone (references 6 and 7) 
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TABLE II - EMPIRICAL COEFFICIENTS FOR CALCULATING MAXIMUM FLAME 


VELOCITY FROM HYDROCARBON STRUCTURE (REFERENCE 26) 



Type C-H bond 


Coefficient 

Methane 

k a 

35. 2X10 " 19 

Primary 

k b 

42.5 

Secondary 

Kc 

47.5 

Tertiary 

Kd 

45.4 

Cyclohexyl 

«g 

50.5 

AiLkene 

k e 

80.7 

Aromatic 

k h 

84.3 

Alkyne 

k f 

223.9 


TABLE III - BOND DISSOCIATION ENERGIES FOR C 4 AND Cg HYDROCARBONS 
AT CONCENTRATION FOR MAXIMUM FLAME VELOCITY 


Hydrocarbon 

Equiva- 
lence 
ratio 
for max- 
imum U, 

^max 

Deviation 

^max 

from average 
%iax 

(percent) 

Dissocia- 
tion energy 
of hydro- 
carbon 
/ kcal \ 

Dissocia- 
tion energy 
of mixture 

^max 
/ kcal \ 

\ liter/ 

Devia- 

tion 

from 

average 

(percent) 

\g-moleJ 

(a) 

Butane 

1.09 

5-5 

1250 

1.93 

2.0 

Hexane 

1.16 

0.9 

1791 

2.00 

1.5 

2 -Me thy lpr opane 

1.11 

3.6 

1257 

1.94 

1.5 

2,2 -Dimethyl- 






butane 

1.12 

2.5 

1821 

1.97 

0.0 

Butene -1 

1.17 

1.7 

1133 

1.98 

0.5 

Hexene -1 

1.16 

• 0.9 

1680 

2.02 

2.5 

Isobutene 

1.14 

0.9 

1147 

1.97 

0.0 

2 -Methylpentene - 1 

1.19 

3.4 

1603 

1.99 

1.0 

Butyne-1 

1.17 

1.7 

1016 

1.95 

1.0 

Hexyne-1 

1.21 

5.0 

1564 

1.97 

0.0 

4 -Methylpentyne 

1.18 

2.5 

1574 

1.95 

1.0 

Benzene 

1.34 

14.2 

1305 

1.96 

0.5 

Average 

1.15 

3.29 

— 

1.96 

0.9 


Single-bond, energies from reference 27. Multiple-bond energies from 
reference 28: 151.2 kcal/g-mole for C=C in olefins and cyclic compounds; 

198.5 kcal/g-mole for C=C bond. 
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.6 .8 1.0 1.2 
Equivalence ratio 

(a) Complete experimental range. 

Figure 1. - Predictions of the effect of equiva 
lence ratio on flame velocity. 
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» 6 .8 1.0 1.2 1.4 

Equivalence ratio 


(t) Limited range. 

Figure 1. - Cor, \uded. Predictions of the effect 
of equivale.ee ratio on flame velocity. 



NACA-Langley - 1-2-53 - 500 



Figure 2. - Empirical 



for equivalence ratio data. 
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